The fate of the a- vacuum 



Hael ColhnsQ R. Holman[[| and Matthew R. Martiri 
Department of Physics, Carnegie Mellon University, Pittsburgh PA 15213 
(Dated: February 1, 2008) 

de Sitter space-time has a one complex parameter family of invariant vacua for the theory of a 
free, massive scalar field. For most of these vacua, in an interacting scalar theory the one loop 
corrections diverge linearly for large values of the loop momentum. These divergences are not of a 
form that can be removed by a de Sitter invariant counterterm, except in the case of the Euclidean, 
or Bunch-Davies, vacuum. 
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I. INTRODUCTION 

The importance of understanding quantum field the- 
ory in de Sitter space, the space-time associated with 
a positive cosmological constant, has been heightened 
by recent observations of both the early and late uni- 
verse. The dramatic results of Wilkinson Microwave 
Anisotropy Probe Q have provided further strong evi- 
dence that the universe underwent a rapid inflationary 
expansion. Both the large-angle anti-correlation in the 
temperature-polarization cross-power spectrum and the 
nearly flat spectral index are consistent with the predic- 
tions of inflation. More surprisingly, the dimming of the 
type la supernovae seen by the Supernova Search Team 
and the Supernova Cosmology Project Q, combined 
with other observations, is yielding a new standard pic- 
ture for the contents of the universe, the largest com- 
ponent of which is a dark energy whose properties are 
consistent with a positive cosmological constant. 

A striking difference between de Sitter and flat space- 
time is the richer vacuum structure of the former. For 
a free scalar field in a Minkowski space, there exists an 
essentially unique Poincare invariant vacuum state. In 
contrast, for a de Sitter background, Mottola 3 and 
Allen p| discovered an infinite family of vacua for the 
quantum theory of a free massive scalar field that are 
invariant under the isometries of de Sitter space. These 
vacua can be parameterized by a single complex num- 
ber, a, and are usually called the a-vacua. Most of these 
a- vacua have a host of peculiar features, such as a mix- 
ture of positive and negative frequency modes at short- 
distances and a non-thermal behavior that violates the 
principle of detailed balance. Only one of these states, 
the Euclidean or Bunch-Davies vacuum, behaves ther- 
mally when viewed by an Unruh detector Q and reduces 
to the Minkowski vacuum as we take the cosmological 
constant to zero. The assumption that the universe was 
at least approximately in the Euclidean vacuum underlies 
the successful predictions of inflation for the calculation 
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of the density fluctuations which produced the temper- 
ature anisotropies in the cosmic microwave background 
radiation. 

Despite their unappealing features, the a-vacua are 
perfectly valid vacua for a free scalar field. If they 
cannot be shown to be unphysical, then their existence 
would undercut some of the robustness of the inflation- 
ary paradigm — we would need to explain how the epoch 
prior to inflation managed to place the universe in the Eu- 
clidean vacuum rather than one of the other infinite fam- 
ily of a-vacua. For example, the regularization needed 
by the energy-momentum tensor even for the free theory 
in the a-vacuum is not generally compatible with that 
needed after inflation 

A complication in formulating quantum field theory in 
de Sitter space is its lack of a well-defined S'-matrix. In 
an interacting theory we have two sources of time depen- 
dence for matrix elements — one induced by any inherent 
time dependence of the background geometry and an- 
other introduced by the interactions. In such a system, 
it is therefore appropriate only to ask time dependent 
questions — to study how a matrix element evolves from 
a given initial state. 

Schwinger [jj and Keldysh [T^j developed a formalism 
to solve for this finite time evolution. In their approach, 
we specify the state of the system at an initial time and 
then evolve to a finite time later. Here, both the 'in' and 
'out' states correspond to the same state and are evolved 
together when we evaluate the expectation value of an 
operator — in effect this formalism evaluates matrix ele- 
ments between two 'in' states. The Schwinger-Keldysh 
formalism is thus ideally suited for studying the behav- 
ior of the a-vacua in the presence of interactions. We 
place the system initially in an a-vacuum and then study 
whether a sensible evolution results. Since the quantum 
field theory only is evolved over a finite interval, our re- 
sults are relevant not only for the more formal question 
of the a-vacua in an eternal de Sitter background but 
also for the phenomenological problem of a finite epoch 
of inflation. 

The methods established here can also be applied to 
any initial state, such as the 'truncated a-vacua' of 
In these vacua, the short-distance behavior of the a- 
vacua is modified either in accord with some specific the- 
ory, such as the stringy uncertainty relation of jl2lll3| . or 
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simply by truncating the a mode functions above some 
energy scale to reflect our ignorance of the new physics 
[Til fiH ITM lift Il7j . We address the formal case of a pur e 
a- vacuum and shall study the truncated case later in [18j . 

In this article, we show that an interacting scalar field 
theory in a general a-vacua contains linear divergences 
which cannot be removed with a de Sitter invariant renor- 
malization prescription. These divergences appear in the 
one loop corrections and are present for arbitrarily weak 
interactions. The specific example we study is the expec- 
tation value of the number of Euclidean particles in an 
a-vacuum. The divergences appear in the high momen- 
tum region of the loop integral. We show that they only 
vanish for the Euclidean vacuum, which is completely 
rcnormalizable . 

The subject of the a-vacua for an interacting the- 
ory has also been recently investigated in related work 
[l?l l2fl | . Both of these works essentially studied the cor- 
rections to the two-point correlation function obtained 
between an 'in' a state and an 'out' state given by the 
a state at a later time. Banks and Mannelli |l9j found 
that the interacting theory in the a- vacuum required non- 
local counterterms while Einhorn and Larsen |20| found 
pinched singularities in the loop corrections. These fea- 
tures provided highly suggestive evidence that the a- 
vacua are pathological in the presence of interactions. 
Some attempts to modify the theory to avoid these prob- 
lems appear in 0, ^ . 

We begin with a review of the de Sitter invariant vacua 
for a free scalar field in Sec. [H] This section also shows 
the form of the Wightman functions in conformally flat 
coordinates. Section II I II derives the expectation value 
of an operator in an interacting theory based on the 
Schwingcr-Kcldysh formalism. In Sec. IIVI we calculate 
the change in the number of Euclidean particles in an 
a-vacuum due to a cubic interaction and show that in 
the presence of this interaction, the expectation value 
is renormalizable for the Euclidean vacuum while an 
unrenormalizable divergence appears for the a-vacuum. 
Section Ivl explores the origin of these divergences in the 
a-vacua in a more general setting. We derive the neces- 
sary conditions for these divergences to arise and show 
how they can appear in a general interacting scalar field 
theory. Section IVII summarizes our results and suggests 
future applications for this formalism. 



II. GREEN'S FUNCTIONS 

In this section we review the rich vacuum structure of 
a free scalar field in de Sitter space |2^|. We derive the 
form of the Wightman function and eventually the Feyn- 
man propagator in conformally flat coordinates. These 
Green's functions will be used later for studying the in- 
teracting theory. 

The most straightforward method for demonstrating 
the existence of a family of de Sitter invariant vacua is 
to evaluate the two-point Wightman function for a free 



massive scalar field in an a-vacuum. For this purpose 
it is useful to use a coordinate system that covers the 
entire space-time. Such coordinates are not, however, 
those best suited for more explicit calculations. There- 
fore, throughout this article we shall study de Sitter space 
using conformally flat coordinates, 



ds 2 = 



di] 2 — dx 2 
H 2 V 2 



(2.1) 



with r\ £ [—00, 0] which cover half of de Sitter space [24j. 
The other half of the space is covered by a set of coordi- 
nates with r\ — ► —77. These coordinates are simply related 
to the standard coordinates used in inflation, 



ds 2 = dt 2 



dx 2 



(2.2) 



through r) = —H 1 e Ht . H is the Hubble constant and 
is related to the cosmological constant by A = 6H 2 . 



A. The Euclidean vacuum 

To an observer capable only of probing length scales 
on which the curvature of de Sitter space is not appar- 
ent, the space-time appears approximately flat. For the 
high energy modes then, this observer can apply the same 
prescription for defining positive and negative frequency 
modes as in Minkowski space. The vacuum state anni- 
hilated by the operators a~ associated with these modes 
corresponds to the Euclidean vacuum. 

The Euclidean vacuum possesses many desirable prop- 
erties in addition to matching with the Minkowski vac- 
uum at short distances or as H — > 0. It corresponds 
to the unique state whose Wightman function is ana- 
lytic when continued to the lower half of the Euclidean 
sphere. Moreover, an Unruh detector placed in the Eu- 
clidean vacuum satisfies the principle of detailed balance 
as though it were immersed in a thermal system at the 
de Sitter temperature, T^s — H/2tt [23j. 

If we denote the Euclidean vacuum by \E), the Eu- 
clidean Wightman function for a free massive scalar field 
$(x) is defined by 



G E {x,x') = (E\$(x)Q(x')\E). 



(2.3) 



Since the metric is spatially flat, the expansion of the 
scalar field <&(x) in creation and annihilation operators 
a-\a£ for the vacuum state \E) is 

k k ' ' 



$(r?,f) 



d 3 k 



(2tt) 3 

With the commutator normalized to be 



-ik-x 

k 



(2.4) 



(2.5) 



the Euclidean Wightman function in position space is 



G E (x, x') 



d 3 k 
(2tt) 3 



Jk-(x- 



3 



d 3 k 

(27T) 5 



e^^Gf (»,,»/), (2.6) 



where the momentum representation the Wightman func- 
tion is 



(2.7) 



Note that the mode functions only depend on the mag- 
nitude of the spatial momentum, k = |fc|. 

A free massive scalar field satisfies the Klein-Gordon 
equation, 

[V 2 + m 2 ] $(x) = 0, (2.8) 

so that the mode functions solve the differential equation, 

[r) 2 d 2 - 2r)d v + ifk 2 + m 2 H- 2 ] Ug(rj) = 0. (2.9) 

Note that m 2 here represents the effective mass of the 
theory which includes any contribution from coupling the 
field to the curvature, <E> 2 i?. In de Sitter space-time the 
curvature R is constant so this coupling is of the same 
form as a mass term. The solutions to Eq. (|2.9|l are linear 
combinations of Bessel functions, 



Uk(v) = c fe V 3/2 Mkrj) + d k 7f/ 2 Y v {kr)), (2.10) 



with 



\-m 2 H~ 2 . 



(2.11) 



We shall assume hereafter that v is real. 

The general form for the mode functions is applicable 
to both the Euclidean vacuum and the a-vacuum. What 
distinguishes the former is that as the Hubble constant is 
taken to vanish, H — > 0, so that de Sitter space becomes 
flat, we should recover only the positive frequency mode 
functions, e~ %kt . In the small H limit, 



krj 



ke~ m 



- + kt + 0(H), 



(2.12) 



the leading time dependence of the modes is U®{rj) oc 



ikl 



when 



-iNu. 



(2.13) 



Up to the normalization factor, Nk, the Euclidean mode 
functions are given by 

Uk(v) = N k rf' 2 \J v (krj) - iY v {krj)] 

= N k r)V 2 HW(kr,). (2.14) 

H{i {kr{) represents a Hankel function. We shall now 
choose the units such that H = 1. 

The normalization is fixed by the canonical equal-time 
commutation relation 



[n(r/,x),$(ry,x')] = -i5 3 {x-x') 



(2.15) 



where the conjugate momentum is 



1 



U{rj,x) = —dn$(ri,x). 



(2.16) 



The equal-time commutation relation requires that the 
modes satisfy a Wronskian condition of the form 

U? d v Ut - d n U k E Ut = i V 2 , (2.17) 

which determines the normalization of the modes to be 



N k 



(2.18) 



Therefore, the Euclidean mode functions are given by 



U k E (v) = ^V 3/2 Hl 2) (kr 1 ). 



(2.19) 



While the de Sitter invariance of the Wightman func- 
tion is not manifest from Eq. H2.19fl . it is possible to 
write Ge(x,x') as a function of the de Sitter invariant 
distance between its arguments 23] . In conformally flat 
coordinates, this invariant distance between x — (77, x) 
and x' = (77', x') is 



Z(x, x) = 



rf + r)' 2 — \x — x'\ 
2rjrj' 



•/|2 



(2.20) 



Although we shall state most of our results in terms of 
the mode functions for a general mass, it will be conve- 
nient to show the results for a particular case in which 
the mode functions simplify substantially. When v = i, 
the Hankel function in Eq. 1)2. 19|) is proportional to an 
exponential, 



/2 



rje 



-ikrj 



(2.21) 



This case corresponds to a massless, conformally coupled 
scalar field for which the effective mass is m 2 — 2. The 
Euclidean Wightman function is then 



G E (x, x 1 ) 



16tt 3 



d A k 
~k~ 



-ik(t]—rj ) ik-(x—x ) 



(2.22) 



and is finite provided we choose the appropriate ie pre- 
scription, 



G E (x,x') = - 



1 



1 



87T 2 Z — 1 — ie sgn(?y — ?/) 



(2.23) 



Here the appearance of the invariant distance Z(x, x') 
establishes the de Sitter invariance of the vacuum. 



B. The a- vacua 

The choice of the short distance behavior of the mode 
functions which determined the relative contributions of 
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the two independent solutions to the Klein-Gordon equa- 
tion is not the unique choice which leads to a de Sitter 
invariant Wightman function. Mottola Q and Allen [f| 
observed that the vacuum state |a) annihilated by a Bo- 
golubov transformation of the Euclidean operators, 



a k 



4 



a -l 



(2.24) 



also yields a de Sitter invariant Wightman function, 

G a {x,x') = (a\${x)§{x')\a). (2.25) 



Here, Re a < and the normalization 



N a = 1-e 



-1/2 



(2.26) 



is chosen to preserve the normalization of the commuta- 
tion relation in the a- vacua, analogous to Eq. I|2.5|l . Note 
that the Euclidean vacuum is itself among the a-vacua 
being obtained when a — > — oo. 

In proving that G a (x,x') only depends on Z(x,x') it 
is useful to use a coordinatization that covers the entire 
de Sitter space-time. In such global coordinates, both a 
point x and its antipode x a occur in the same coordinate 
system 23] . It is then possible to chose Euclidean mode 
functions 0^(x) such that cj)^*(x) — ^(xa) so that the 
Bogolubov transformation of Eq. 12.24fl gives 



<%(x)=N a [ct>%(x) + e a <t>Z(x A j] 



(2.27) 



Here n labels the elements of a general basis of mode 
functions. In this form, the de Sitter invariance of the Eu- 
clidean Wightman function and the fact that Z(x, x' A ) = 
— Z(x,x') together imply that the a- Wightman function 
only depends on the de Sitter invariant distance between 
x and x 1 [4J, |5j . While it is helpful to use a coordinate 
system which contains the antipode of every point to es- 
tablish this invariance, for explicit calculations it is not 
necessary to use global coordinates. Equation (|2.24|l re- 
lates the mode functions of the a- vacuum to the mode 
functions of the Euclidean vacuum and their complex 
conjugates — we do not need to transform the conjugated 
mode function into a function of the antipode once we 
have established that G Q (a:,a;') is invariant. 

From the Euclidean mode functions of Eq. I|2.19|l we 
can now construct the mode functions for the a-vacua. 
Expanding the scalar field in terms of a creation and 
annihilation operators, 



$(77, £) 



d 3 k 
(2^)3 



UZ(r,y k -*a1 + Ur(ri)e 



-ik-x ^af 



and using Eq. (|2.24|) yields 

U£(r,) = N a [Uj?(r,) + e a Ut(ri)] 



k 

(2.28) 
(2.29) 



since the Uj?(rj) only depend on the magnitude of k. Thus 
the a- vacuum modes are 

U?(V) = N a ^/ 2 iH^ikrj) + e"HW(kr,)] . (2.30) 



Inserting the a-mode expansion into Eq. (|2.25() yields 



G a {x,x') = 



d 3 k 
(2tt) 3 



iki*-Z') mr])u «* {r] , )t (2 _ 31) 



Again, the spatial flatness makes it natural to use a mo- 
mentum representation 



G%( V ,r 1 ') = U%(r ) )UrW)- 



(2.32) 



The additional complexity of the mode functions in the 
a-vacuum means that it is particularly helpful to have a 
case in which these functions simplify. For a massless, 
conformally coupled scalar field, 



ikrj 



and the Wightman function becomes 
N l f 1 

G a (x,x) = — — jr< — ; 

87H L Z — 1 — ie sgn 



e a e ik71 ] 



(2.33) 



(v - v') 



e 



Z — 1 + ie sgn(f] — ?/) 

e a e a * 



Z+l 



z + 1 



(2.34) 



As in the Euclidean case, the de Sitter invariance is man- 
ifest in the above expression. 



C. Propagation 

To study the propagation of signals in a de Sitter back- 
ground, define the Feynman propagator, 



iG(x,x') = (a\T($(x)<f>(x'))\a), 



(2.35) 



so that it satisfies the Klein-Gordon equation with a point 
source, 



[V^+m 2 ] G(x,x') = 



S 4 (x - x 1 ) 



(2.36) 



The propagator can only depend on the difference be- 
tween the spatial positions of the points so that its 
Fourier transform is 

G k (v,v') = I d 3 xe^< s ~^G( V ,x; V \x') (2.37) 



where Gk(rj,r)') is the solution to 

[jfdl - 2 V d v + ifk 2 + m 2 ] G k ( V , rf) 

= r] 2 r]' 2 5(r]-r)'). (2.38) 

The solution to this equation which satisfies the correct 
boundary conditions at 77' = rj has the form 

G k (v, V') = G>( V , t/) 6(77 - r/) + G<( V , rf) 6(77' - 7,). 

(2.39) 
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Here G k (77,77') an d (77, 77') are essentially the 
two-point Wightman functions calculated earlier in 
Eq. 

G>(v,v') = iG%{v,r/)=iUgWk*tf) 
G<(r,,r,') = iG%{r 1 ',ri)=iUr{r,)U%{r ] '). (2.40) 

Although propagation for a free field theory in the a- 
vacua contains some peculiar features, it is not otherwise 
ill-defined. The pathological features of quantum field 
theory in an a- vacuum only appear in an interacting field 
theory. The form of the a- Wightman function already 
suggests that the interacting theory could be ill-defined, 
since the various terms in Eq. I|2.34|l contain different 
ie prescriptions. This property implies that in a stan- 
dard approach to calculating the one loop corrections to 
the propagator, among the products of the propagators 
participating in the loop appear products of poles with 
the opposite ie prescription — pinched singularities |20| . 
For example, in the one-loop correction to the propa- 
gator appears a product of the Green's functions given 
in Eq. I|2.34|) . However, these pinched singularities do 
not by themselves prove whether the a-vacuum is itself 
pathological or whether the standard methods for study- 
ing the quantized theory are inappropriate for a time- 
evolving background such as de Sitter space. 



III. THE SCHWINGER-KELDYSH 
FORMALISM 

A significant difference between de Sitter space, in 
most coordinatizations, and flat space is the explicit 
time-dependence of the metric. Unlike a flat space-time 
where the generator of time-translations is a Killing vec- 
tor globally, de Sitter space-time has no such global 
time-like Killing vector. Moreover, in a particular co- 
ordinate system — such as inflationary coordinates — the 
time derivative may not even generate an isomctry lo- 
cally. These properties suggest that rather than attempt 
to define an 5*-matrix between 'in' and 'out' states de- 
fined at different times, we should apply a quantization 
procedure that evolves an entire matrix element over a 
finite interval. It is also useful to be able to evolve a given 
state forward from a specified initial time 770, rather than 
to use a state in the asymptotic past. We can always 
take 770 - * —00. 

An additional advantage of solving the evolution over 
finite intervals is that such an approach more immedi- 
ately determines whether the a-vacuum is applicable for 
inflation, which does not require a de Sitter space-time 
eternally, but only over a sufficient interval to generate 
the number of e-foldings needed to explain the flatness 
and the homogeneity of the universe. If the interacting 
a-vacuum shows its pathology even over a finite interval, 
then we can exclude the possibility that the universe was 
in a pure a-state during any epoch of inflation, regardless 
of the prior history of the universe. 



The closed time contour formalism developed by 
Schwinger 0, Keldysh an d Mahanthappa [2^ al- 
lows us to study the evolution of a quantum field theory 
over a finite interval after specifying the state at an ini- 
tial surface. We review here their approach which leads 
to an expression for perturbatively evaluating the matrix 
element of an operator, which is given at the end of this 
section in Eq. l|3.I5fl . 

In the interaction picture, the evolution of operators is 
given by the free Hamiltonian, Hq, while the evolution 
of states is given by the interactions, Hj. If we let {| V E')} 
denote a general basis of states for the theory, then the 
behavior of the system is completely described by the 
density matrix, p{rf) — ^, p*,*' 1^) (^'l- Thus, as the 
density matrix is constructed from the states, it satisfies 
a Schrodinger equation of the form 

i—p(r l )=[H I ,p(r 1 )]. (3.1) 

The advantage of the interaction picture is that fields 
evolve using the free Hamiltonian, 

-i^( V ,x) = [fT , *(»?,*)]• (3.2) 

The time evolution of the field is precisely that given in 
the previous section since here the mode functions still 
are solutions to the free Klein-Gordon equation. 

To study the evolution introduced by the interactions, 
it is convenient to include a 'turning on' function in the 
interacting part of the Hamiltonian, 

H = H +oj(r ) -r lo )H I ; (3.3) 

here uj(r~i — 770) vanishes when 77 < 770 and becomes one 
when 77 is sufficiently large compared with 770. Later we 
shall let this function be a O step function. We shall 
often not write this function explicitly, absorbing it into 
Hj. Thus the state does not evolve before 770: p(rj) = 
p(Vo) = Po for 77 < 770. 

Once we have specified the state at a particular time, 
^(770), then the Eq. I|3.I[) allows us to determine the state 
at all subsequent times. To study the vacuum structure 
of de Sitter space, the initial state will correspond to 
an a-vacuum. If we introduce a unitary, time-evolution 
operator £7/(77,77') that evolves the state, 

P(jl) = U I {r],r]o)p{vo)Ui 1 (ri,rio), (3.4) 

then from Eq. I|3.1|l U (77, 770) obeys 

d 

ig-Ui(r], 770) = ii/£//(77, 770) (3.5) 

with [/f(?7o, 770) = I. The formal solution to this equation 
is given by Dyson's equation in terms of the time-ordered 
exponential 

^(77,770)= Te- iJ ™ (3.6) 
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FIG. 1: The contour used to evaluate the evolution of op- 
erators over a finite time interval. The initial state is an 
eigenstate of the Hamiltonian until 7,0 at which time the in- 
teractions are turned on. We double the field content so that 
separate copies of the fields are used for the upper and lower 
parts of the contour. 



The evolution of the expectation value of an operator 
in this time-dependent background is given by 



(0)(v) = 



Tr[p(r?)Q] 
Tr[p(r,)] 

Ti[p Uf\ v , r^OUiiju 7,0)] 
Tr[po] 



(3.7) 



Since the state p(r)o) does not evolve before the inter- 
actions are turned on, we can insert the identity in the 
form Ui(r]o,r]p)Ui(rip, 770) with r\ v < 770 and commute one 
of these evolution operators with po to obtain 



Tr [poUiijfr, if)OUi{rj, rjpj] 
TrU] 



(3.8) 



Inserting another factor of the identity, 
Ui(v P ,Vf)Ui(f]f,Vp), witn Vf > V yields 

, . / x = Tr[p o t/ J (??p,7? / )t/j(?? / ,^0L r /(??,^)] ^ ^ 
Tr[po] 

and finally we let rj p — > — 00 and rjf — ► 0, which rep- 
resent the infinite past and infinite future in conformal 
coordinates, so that 



(o)(v) = 



TrjUij-^, 0)Ui{0, v)OUi(v, -00) Po] 
Tr[Ui{-oc,0)Ui(0,-oc)p ] 



(3.10) 



Reading the operators from right to left in the numerator 
of this equation, po sets the initial state of the system 
which is then evolved along a time contour from —00 
to with an operator inserted at 7,; the final operator 
evolves back from to —00. The closed time contour 
which results is depicted in Fig. ^ 

To evaluate Eq. (|3.1U|I it is useful to group the evo- 
lution operators into a single time-ordered exponen- 
tial. This is accomplished by formally doubling the field 
content of the theory, with a set of '+' fields on the 
increasing-time contour and a set of ' — ' fields on the 
decreasing-time contour. The arrows on the contour in- 
dicate time ordering of events so that events on the — 
contour always occur after those on the + contour. We 
can group the effects of both parts of the contour together 
by writing the interacting part of the action appearing in 
Dyson's equation, as 



Si = - 



d?, #/($+) - 



d-nH^®-). (3.11) 



Since the two terms differ only in the direction of the 
integral over the conformal time, we can write the action 
as a single Lagrange density, 

Si = - f dq [Hi(<P+) - Hi($-)] . (3.12) 

J — 00 

The field doubling induced by the closed contour ef- 
fectively doubles the number of vertices we must include 
when studying any process — one set with fields on the 
+ branch and one with fields on the — branch. From 
the Eq. 1)3. 12[) . the latter will have couplings with the 
opposite sign. In evaluating matrix elements, Wick con- 
tractions produce four propagators for the possible con- 
tractions of pairs of the two types of fields, 



(alT^in, f)* ± (j 7 ',i'))|o) 



d 3 k 



lk-(x — X ) S~1±± 



GfW). (3.13) 



The time-ordering of the contractions is determined by 
the direction along the contour, 

C++ (77, 77') = G> (77, 77')6(77- 77') + G< (77, 7,06(7,' -77) 
G^~M) = G> (77, 7,06(77'- 7,) +G< (77, 7,06(77 -77') 
G fc -+(7,,77') = G>(7,,7,') 

G+-(V,V') = G<(t,,7,0, (3.14) 

with the Wightman functions given in Eq. I|2.40|l . 

Assembling the ingredients of the Schwinger-Keldysh 
formalism — the general expression for an operator expec- 
tation value in Eq. (|3.10() . Dyson's equation l|3.6|l and 
Eq. H3.12f> — provides an explicit expression for the evo- 
lution of (0)(rf). If we let the initial density matrix be 
that for a pure a- vacuum, then Eq. (|3.1U|I becomes 



(a\0\a)( V ) 



(a\T{ 



Oie 



}l«) 



(a|T{« 



}l«) 



... ( 3 ' 15 ) 
Here we have absorbed any 'turning on' function in Hi — 

in essence the time integrals begin at 7,0. The time or- 
dering has allowed us to group the time evolution oper- 
ators in Eq. I|3.10|l along the two contours into a single 
operator. This equation for the finite evolution of the 
expectation value of O is the analogue of the standard 
5-matrix expression used in Minkowski space. The virtue 
of the field doubling is that it removes any acausal be- 
havior from the matrix element since the 6-functions in 
the propagators combine to limit the upper end of the 
conformal time integrals to 7, 0, llfjfl . This property will 
become clear in the calculation of a specific example. 



IV. EVOLUTION OF THE NUMBER 
OPERATOR 

We now evaluate the expectation value of the num- 
ber of Euclidean particles in the cv-vacuum using the 
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Schwinger-Keldysh formalism |26(. This number oper- 
ator provides a good measure of whether a particular 
choice for the vacuum state becomes pathological in the 
presence of interactions. From the perspective of the Eu- 
clidean vacuum, the a-vacuum is an excited state. We 
can determine the stability of this state when the inter- 
actions are turned on by following the time evolution of 
the change in the Euclidean number operator evaluated 
in the a-vacuum. 

In the non-interacting theory, we should not encounter 
any infinite Euclidean particle production in the a- 
vacuum. When the interactions are turned on, some 
further particle production will occur in the Euclidean 
vacuum, but the rate per unit volume should be finite 
and any divergences which appear perturbatively must 
be renormalizable. What we shall discover is that in the 
interacting theory, the a-vacuum produces a new class of 
divergences that cannot be removed with the usual set of 
counterterms. 

These divergences in the a-vacuum occur at each time 
in the integrand as we propagate from some initial state 
at 770 to a finite time later so that the theory diverges even 
for an arbitrarily short time after the interactions are 
turned on. This behavior indicates that the a-vacuum is 
inappropriate even for a finite inflationary epoch. 

Although we examine in detail the evolution of the Eu- 
clidean number operator in this section, all of the diver- 
gences we find are generic to the one loop corrections to 
an arbitrary operator. Furthermore, while we consider 
a scalar theory with a cubic interaction since it has a 
simple, non-trivial self-energy correction, similar loop in- 
tegrals occur in any interacting scalar field theory. These 
more general cases are treated in the next section. 

We first show how the free Hamiltonian can produce a 
non-trivial, but finite, time dependence. We then set up 
the calculation for the expectation value of the deriva- 
tive of the number operator, expressing it in terms of 
the scalar field and its conjugate momentum. The cor- 
rections to this expectation value to leading, non-trivial 
order in the coupling are then calculated and we obtain 
general expressions for the one-loop corrections and the 
counterterms. 

In the interaction picture, the Hamiltonian is divided 
into free and interacting parts, 



H = H + e( V -r l0 )H I . 



(4.1) 



Here we have included a O-function so that before 770 the 
system evolves freely; we are always free to take 770 — > 
—00. For conformally flat de Sitter coordinates, the free 
Hamiltonian for a scalar field $(77, x) of mass m is given 
by 



H = / d 3 x 



Vn 2 + ±?r 2 (V$) 2 + ±77- 4 m 2 $ 2 . (4.2) 



An important difference between de Sitter space and flat 
space is that the free Hamiltonian is not diagonal in terms 



of creation and annihilation operators 
d 3 k 



rr _ 1 1 



(2tt) 3 



E E\ 



M E 
a k a k 



9k (v) 



+afa^{r 1 )+afa E U E *{r ] )\ (4-3) 



where 



f E 
Ik 



= d v u£d v u£ 



fc 2 + -jt 



9 II V 

fc 2 + 



>r 



u k E u h E . 



(4.4) 



Note that this property holds also for the a-vacuum, E — > 
a. 

The origin of the off-diagonal terms in the free Hamil- 
tonian lies in the fact that the surfaces of constant con- 
formal time are not orthogonal to the generator of an 
isometry. This effect introduces an additional non-trivial 
source of time evolution which combines with that pro- 
duced by the evolution of the state |a) when the fields 
interact. For a time derivative of a generic operator, we 
formally have 

9,(0(77)) - Tt[{d v p)0 + pd v O] 

= iTr [-[Hi,p]O + p[H ,O]] 

= iTr[p[H,0]} =i{[H,<D]}. (4.5) 



The off-diagonal terms in the free Hamiltonian induce 
an evolution in the number operator even in the free the- 
ory. Since we wish to explore the effect of interactions on 
the a-vacuum, we construct the number operator from 
creation and annihilation operators, at and &£, which 
satisfy 



4 



(4.6) 



at the moment the interactions are turned on. The sub- 
sequent evolution in the interaction picture is given by 
the solution to 



.d_ _ 

d?7 k 



(4.7) 



The solution to this equation can be formally expressed 
as 



(4.8) 



where Uq(j], rjo) is the time evolution operator for the free 
part of the theory, 



U (r),no) 



Te -iJZ dv"H ( v ")_ 



(4.9) 



From the form for the free Hamiltonian, a general solu- 
tion to Eq. (|4.7|l is given by a Bogolubov transformation 
of the time independent creation and annihilation oper- 
ators, 



a k (v) a,: + /3 fc (7?) a f -., 



(4.10) 
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where the coefficients satisfy 

ir] 2 d v a* k = #f a£ + & 
-irj 2 d v /3 k = g®(3 k + fk* a t- 



(4.11) 



The standard normalization of the commutator of the 
transformed creation and annihilation operators also re- 
quires that 



\a k (v)\ 2 -\Mv)\ 2 = l- 



(4.12) 



The general solution to the coefficient equations (|4.11|) 
is of the form 



a 2 



where the constants a\ and a 2 should satisfy 
aia^ — a\a 2 = —i 



(4.13) 



(4.14) 



from Eq. H4.12[l . If we would like the number operator to 
count the number of Euclidean particles at the moment 
we turn on the interactions, 77 = t]q, then 



at 



V 2 o 



a 2 = iUtiVo) 



(4.15) 



The number operator constructed from the trans- 
formed creation and annihilation operators, 4^fc' has the 
correct evolution for an interaction picture operator. Us- 
ing Eq. (|4.5|l . even the free Hamiltonian induces some 
evolution in this number operator, 



Evaluating Eq. (|4.16() in the a-vacuum gives the tree 
level evolution of the number operator in this state, 

i(a\[H ,a)ja s ]\a) = -^VD ab (r])U^(r] a )U^*(r] b ) 

(4.19) 

V is the spatial volume. We can divide the volume from 
both sides to yield the particle production rate per unit 
volume for which this tree contribution is completely fi- 
nite, if non-zero. 

Now consider a cubic interaction with its associated 
counterterms, 

H I= J d 3 xrT 4 [J$ + ±<W$ 2 + ^A$ 3 ] . (4.20) 

The cubic vertex will generally introduce one-loop correc- 
tions to the two-point functions, so we have included a 
mass counterterm, 8m 2 . Since the <I> 3 interaction breaks 
the $ — > — $ symmetry of the free theory, we also expect 
the interaction to generate graphs containing tadpole in- 
sertions which are cancelled with the correct choice for 
J. To the order we shall study, no wave function renor- 
malization is needed. 

The change in the number operator induced by these 
interactions is given by applying Eq. (|3.15|) . 

i{a\T { [H, 4 a % ]e- { C *> W+Wn] } | a) 



(a\T{ 



f !(*-)] 1 



[ffo,4«d = — D ab J d 3 fd 3 ye lfc ^$(?7 a ,f)$(^,y), 

(4.16) 

where D a b(jf) is the differential operator, 

D ab {n) = Pl ( v ) 3 Va d nb -p 2 (r)) (8 Va +d m )+p 3 (r ) ). (4.17) 

Here r) a ^ are only labels to indicate the functions on 
which the derivatives act; in the end, r\ a ^ are set equal 
to 77. The functions Pi{rj) are 



pi(v) = 






Vo 




Pa(»?) = 


-Ad,ui{ m )\ 2 








P3(J?) = 




, 2 ™ 2 
k H 5- 




Vo 





, 2 m 
T] 



(4.18) 



The evolution of the free field is simple in the interaction 
picture, so it is useful to write the creation and annihi- 
lation operators in terms of the field and its conjugate 
momentum, by expanding in the time independent oper- 
ators, Eq. (|4.1U|I . and inverting the operator expansion 
in Eq. ij^Ij) . 

af = i [ d 3 xe~^ [Utn(v,x) - V - 2 d v UtHv,x)} 



k 



(4.22) 



The commutator with the free part of the Hamiltonian 
is then as was given in Eq. H4.16(l while that with the 
interacting part is 



rJ(2 



d 6 x 



2\U^(r l0 )\ 2 U(r ] ,x) - r, 2 d v \U* M\ 2 HV^) 



(4.23) 



rf 2 



--^m 2 I d 3 xd 3 ye lk < s -y~> \U^( Vo )\ 2 [$(77, x)IL(r,, y) + 11(77, f)$(r?, y)] 



d 3 xd 3 y^< s ~^ \\U^( V0 )\ 2 [$ 2 (V, m(v, V) + n(r?, x)$ 2 (77, y)] 



I 




./ 

• x 



= 



.4/ 




FIG. 2: The coefficient J of the linear counterterm is chosen 
to cancel insertions of tadpoles. The dashed line represents a 
line in a general diagram. 



The expectation values of each of these commutators will 
be evaluated in the a- vacuum to order A 2 . 

Before evaluating the expectation values of these com- 
mutators perturbatively, a large class of graphs, those 
containing a tadpole subgraph, are eliminated through 
the proper choice of the coefficient J of the linear coun- 
terterm. To leading order in A, this choice for J is 



J = - 



A 



I6tt 3 



d 3 p \U° (77") I 



(4.24) 



This cancellation is shown diagrammatically in Fig. 
Note that while the loop integral in Eq. (|4.24|) contains an 
apparent time dependence, it is in fact time independent, 



= N 2 TT 2 



ffW(0+e a i#>(0 .(4.25) 



The form of the leading corrections to the expectation 
value of [Hi, atari in the a- vacuum are simpler since the 

k ™ 

commutator is already itself of order A. The order A 2 
corrections to this commutator are shown diagrammat- 
ically in Fig. |3 The other corrections from the cubic 
interaction at this order contain tadpole subgraphs and 
are cancelled when Eq. (|4.24|l is satisfied. 

One of the subtleties in evaluating the evolution of the 
number operator is that it contains time derivatives. For 
example, the Sm 2 term of Fig. [31 contains a term of the 
form 

(a\T ($(77, x)H(r], y) + 11(77, 2)) H (4-26) 

which can produce Schwinger terms if the time ordering 
of the operators is not treated carefully. A method for 
avoiding such terms, following |26j |. is to write the canon- 
ical momentum as 



11(77, x) = jim -f^dn^H 1 ,x) 

77" -H7 77 



(4.27) 



FIG. 3: After choosing the linear term to cancel graphs con- 
taining a tadpole, only these terms contribute at order A 2 
to the expectation value of [i/r,atag] in the a- vacuum. The 
first represents a self-energy graph while the second is from 
the mass counterterm. 

and then to place the fields on the appropriate contours 
so that the time ordering naturally given along the con- 
tour preserves the correct ordering of the operators in 
Eq. KZty . 

lim -4^»(a|T($-(77,f)$+(77",y) + 

$-(?7",f)$ + (?7,f))|a). (4.28) 

With this prescription, the mass counterterm in Fig. 
contributes 

Cj = -^v{\u?to>)\*a,\mri)\ a 

-^d^Mfmirj)] 2 ) (4.29) 

at lowest order. 

The only non-trivial effect of the interactions on ex- 
pectation value of the commutator with the interaction 
Hamiltonian arises from the self-energy graph, 

* = -$Zs\u k E M\ 2 



x 



^Im[cW(^r(^(W)] 



A 2 V dr,\UF( m )\< 



n 4 8tt 3 



f 1 dn' 

x / 4lm[C/ fc a (77)t/^*(77')^(77 , 77')] (4.30) 
where the loop integral is given by 

LUv,ri')= [ <??U?(Ti)U2r(Tf)UZ_ k (T,)U£ h (Tf). 



(4.31) 
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A 



Co 




5m 2 



FIG. 4: The order A 2 corrections to the expectation value of 
[Ho, a? a,j-] in the a-vacuum. Again, the first represents a 
self-energy graph while the second is from the mass countert- 
erm. 



In addition to the corrections in Eq. H4.29|) and 
Eq. (|4.30l) . the fact that the free Hamiltonian is not the 
conserved quantity associated with a time-like genera- 
tor of an isometry of de Sitter space means that the ex- 
pectation value of [Ho,ala£\ also contributes at order 
A 2 . The linear counterterm also cancels the graphs con- 
taining a tadpole subgraph and any vacuum to vacuum 
disconnected graphs are removed by the denominator of 
Eq. H4.21fl . The only remaining corrections then are those 
shown in Fig. These graphs contribute 

A 2 V 



A ° = — iTl^W 

?7 4 27T J 



J no J -no 7 fe 

xIm[UZ(? lb )Ur(m)L%(VuV2)} (4.32) 



and 

2Sm 2 

Co = -^VD ab ( V ) 

"dm 



Im [U£ fa ) Ur ( m ) U% ( Vb )Ur ( Vl )] (4.33) 



to the expectation value of the derivative of the number 
operator, N E 

a,k 

An important feature to note is that the mass coun- 
terterm insertions, in Eq. 1)4. 29|) and Eq. I|4.33|l . do not 
vanish in the Euclidean limit; they are needed to cancel 
a logarithmic divergence in the corresponding self-energy 
diagrams. 



A. Renormalizing the Logarithmic Divergence 

The one loop corrections in the generic a-vacuum and 
in the Euclidean limit differ in their divergence structure. 
Both cases contain a logarithmic divergence which can be 
regularized and cancelled by the appropriate choice for 
the mass counterterm, Sm 2 . What distinguishes the gen- 
eral a-vacuum from the Euclidean vacuum is the appear- 
ance of an additional class of terms that diverge linearly 
in the loop momentum. This divergence cannot be renor- 
malized. In this subsection, we summarize the renormal- 
ization of the logarithmic divergence for the case of the 



massless conformally coupled scalar field in the Euclidean 
vacuum. The detailed calculation for the a-vacuum is left 
for the Appendix lAl 

Formally, the evolution of the Euclidean number oper- 
ator to order A 2 in the Euclidean vacuum is given by 



-^VDMU£(ria)Ut(Vb) 

lim [Ai + d +Ao+C Q ] (4.34) 



Q — i — DC 



The case of a massless, conformally coupled scalar field 
is most readily analyzed since the mode functions have a 
simple form given in Eq. (|2.21() . The tree contribution in 
this case is given by 



V 



2fcV 



Vo Vo no if. 



(4.35) 



which is finite for 77 £ [770, 0). The loop integral for v = i, 
which occurs in Ai, is 



= iir {rjrff ik(ri-ri') 



2 T) - 77' 



(4.36) 



so that 
Ai = 



A 2 


V 




1 


r h 


k 2 


32tt 2 


y 


VoV_ 


J no V-rf 



A 2 V v 2 



dri 



k 32ir 2 77° j i] — T] 1 
The corresponding counterterm is 
Sm 2 V 



— cos2fc(77 — 77') 



7 sin 2^(77-77'). (4.37) 



Ci = 



k 2 2 



77 s 



(4.38) 



The second integral in Eq. (|4.37|l is completely finite at 
all times 77 < 0. The first integral, however, contains a 
logarithmic divergence at the upper end of the drj 1 inte- 
gration. This term can be regularized as described in the 
appendix^ The pole as the regularization is removed, 
e — * 0, is given by 



Ar = 



A 2 V 
k 2 32^ 



4 
77 s 



1 



— h finite, 



and is cancelled by choosing 



Sm 



1 A 2 



e 16tt 2 



(4.39) 



(4.40) 



in Eq. (|4.38|) . The analogous logarithmic divergence in 
Ao is cancelled by the counterterm Cq. 

The self-energy diagrams in the a-vacuum also con- 
tain a logarithmic divergence which can be removed by 
a suitable choice for the mass counterterm, 



Sm 



1 A 2 



e 16tt 2 



1 



N 



(4-41) 



which reduces to that for the Euclidean case in Eq. 14.40fl . 
The origin and regularization of this divergence is dis- 
cussed more fully in Appendix lAl 
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B. The Linear Divergence of the a- vacuum 

Including the order A 2 corrections, the derivative of the 
number of Euclidean particles in the a vacuum is again 
given by the sum of the contributions shown in Fig. Inland 
Fig. 0] as well as the tree level contribution of Eq. I|4.19|l . 



i 



V 

-A! 



d + Aa+Co. (4.42) 



Each self-energy graph contains a loop integral, 
Eq. (|4.31|l . Unlike the Euclidean case which is completely 
finite once we have established an ie prescription, the 
loop integral over a-mode functions diverges linearly in 
the spatial momentum. Introducing a bound A to re- 
move the large momenta in the loop, J °° p 2 dp J dfl2 — > 

p 2 dp J dfl 2 , assuming A > |fc|, the divergent part of 
Eq. l|OT|l is 



k 

sin k(rji 



{V1V2) 



(4.43) 
sinfc(?7i + 772) 



Vi - m 



m + 7] 2 



finite. 



The appearance of the factor e a+a shows why such a 
divergent term does not arise in the Euclidean limit. 

Unlike the logarithmic divergence, this divergence can- 
not be removed by a momentum independent value for 
Sm 2 . For example, the divergent piece of the self-energy 
graph in Fig. is 



A, 



(4.44) 



Hdrf 



sin k(r) — rf) sinfc(?7 + rj') 



+e 



a-\-a 



N- 



r\ — rj' rj + rf 

x 2 vA d^jno)] 2 



71 dp' 



a krj 2 Aw 2 
|lm{^)Lf *(?/)} 
sin k(rj — rf) sinfc(r; 



T] — Tj 



which cannot be cancelled by a momentum-independent 
choice for Sm 2 in Eq. (|4.29p . 

To observe this incompatibility of the momentum de- 
pendence of the divergences in the a-vacua loop correc- 
tions and the available counterterms, at least in the case 
of a massless, conformally coupled scalar field, it is suf- 
ficient to compare the limiting behavior in the external 
momentum of the loop corrections and the correspond- 
ing counterterms. In the k = — ► limit the divergent 
parts of both Ai and Ai scale as fc _1 A. In contrast, the 



leading fc-dependence of the counterterms scale as k~ 2 , 
at least when a is real, which is required in any case for 
a CVT invariant theory (2^. Thus, no choice of Sm 2 , 
which does not depend on k, is possible such that Ci and 
Co cancel the divergence in Ai and Ao as A — > 00. 

Note that the linear divergence is not present in the 
opposite quantity — the expectation value of the number 
of a-particles in the Euclidean vacuum. This quantity is 
given by expressions similar to those above except with 
the labels interchanged for the mode functions, <-» 
Ug , and some of the a-dependent coefficients are slightly 
altered in Eqs. (|4.29|l - (|4.33|) . The crucial difference is 
that the presence of the Euclidean states in Eq. (|3.15|) 
leads to Euclidean propagators so that in particular the 
loop integral is over Euclidean modes, as in Eq. (??), for 
which no linear divergence occurs. 



V. DISCUSSION 

The linear divergence that arises from the one loop 
corrections in the a- vacuum is a UV effect. At arbi- 
trarily short distances there exists an interference of the 
positive and negative frequency modes which cancels the 
rapidly oscillating phases among some of the terms within 
the loop integral. Without such a cancellation, these 
phases could damp these high-momentum contributions 
through an appropriate ie prescription. This interference 
of phases is a specific feature of the propagator in the 
a-vacuum and does not occur in the Euclidean case. In 
this section we shall discuss the origin of this divergence 
and determine the necessary conditions for it to arise. 

Consider a loop containing n vertices connected by n 
internal propagators — those through which the common 
loop momentum flows. Since in de Sitter space it is conve- 
nient to perform a Fourier transform over only the spatial 
coordinates, each vertex has a time, rji for i = 1, . . . ,n, 
associated with it. Eventually, we integrate over all these 
times as they arise from the exponent of the time evo- 
lution operator in Eq. I|3.15fl . To determine whether a 
particular loop can produce a UV divergence, we must 
first count the powers of momentum in the high loop 
momentum region. 

Let Gp_ fc . (77^, 77i + i) represent the Wightman function 
within a loop propagator connecting the i and (i + 1) 
vertices. The loop momentum is p and fcj denotes other 
momenta following through the z-th leg. In the UV limit 
\p — ki\ ^> |77i| -1 , |?7i + i| _1 , the leading behavior of this 
Wightman function is 



(5.1) 



iN 



2 mm+i 

C 2\p-ki 



e —i\p-ki\(rji-rn + i) , e a+a" gi|p-fci|(7Ji-7Ji+i) 



Note that the propagator contains factors of both 
G M,fe.*+i) and G<_ fci (r/ 4 ,?7 2+ i) = G>_ k .(>m + i, rji). 



12 






FIG. 5: Examples of divergent diagrams in theories with quar- 
tic (left) or quintic (right) interactions. Any loop that con- 
tains only two lines will exhibit a UV divergence similar to 
that in Eq. (14. 4311 . This result follows from the structure of 
the propagator in the a-vacuum and not the form of the in- 
teractions. 



For the purpose of the power counting, it is important 
to note that, aside from the phases, in the UV limit 
Gp_ ki (r/i, rji + i) r~j p~ 1 . Thus in integrating over a loop, 



/A n .f 

»=i j 



dp 

pn-2 



(5.2) 



We only encounter a possible UV divergence if n < 3. 
Note that the n = 1 case can be removed by a countert- 
erm since the loop only depends on the loop momentum 
and not on any other momenta in the graph. 

The n — 2 case can produce a linear divergence. 
Since the divergence only depends on the form of the 
propagator and not the form of the interaction, such 
divergences generically occur in any interacting the- 
ory, for example in the processes shown in Fig. 0] and 
Fig. |SJ Superficially, a divergence might seem possible 
even in the Euclidean case if a product of Wightman 
functions, G> (rjx, r)2)G^_ k (r]i, 772), occurs in the loop in- 
tegral. However, the Schwinger-Keldysh formalism is 
constructed to remove such terms and only factors of 
G>(m,m)G> k (r] U V2), and G<(r}i, r/ 2 )G , <_ fe (i7i, r? 2 ) can 
appear in the loop integral in which the p-dependent 
phases do not cancel in the UV. The important difference 
in the a-vacuum is that each Wightman function con- 
tains terms whose p-dependent phases have the opposite 
signs. Thus even after the Schwinger-Keldysh formal- 
ism has been applied, a product of the two loop Green's 
functions, 



G>(r)i,r) 2 )G>_ k (rii,r)2) - 



(5.3) 



a a+a 



N, 



Ap\p — k\ 

1 e ip(m-m) e -i\p-k\(vi-V2) 

_|_ e «p(>Ji +V2 ) e —i\p—te\ (171 +rj 2 ) 
^p-ipim+m) e i\p-k\(rn+r) 2 ) _|_... 

will have some phase cancellation as p — > 00. What ren- 
ders these terms unrenormalizable, however, is that al- 
though the p-dependent phases cancel in the divergent 
terms, they still contain a non-trivial dependence on the 
momentum entering the loop from the rest of the dia- 
gram. 

This analysis of the phase structure of the two- 
propagator loop does not necessarily show that such di- 
vergences cannot be removed from the theory through a 



FIG. 6: Diagrams containing a loop with three propagators 
can be logarithmically divergent in the loop momentum. For 
example, a vertex correction in a "I> 3 theory could generate 
such a divergence. 



suitable rcnormalization prescription. This fact can only 
be established by summing all the contributions to this 
graph and demonstrating that the dependence of the re- 
sulting divergent term on the momenta external to the 
loop is incompatible with a counterterm insertion, as was 
done in the previous section. However, since all loop in- 
tegrals containing only two propagators have essentially 
the same structure, given by Eq. (|4.31|) . we see that the 
a-vacuum cannot be renormalized in any interacting the- 
ory, regardless of the form of the interaction. 

The power counting argument indicates that a loga- 
rithmic divergence can arise for a loop with three legs, 
such as the vertex correction graph shown in Fig. |SJ The 
loop integral contains a product of three terms of the 
form of Eq. (|5.1(l . or its complex conjugate, which gener- 
ally contains terms where the high loop momentum de- 
pendence of the phase cancels. For example, in Fig. 
if the coordinates associated with the three vertices are 
(771, fci), (772,^2) an d (773,— hi — £2), then in the product 
of three propagators occur terms such as 



G>_ fel (r]i,m)Gp +k2 (v2,Vs)Gp (771,773) 

(?7l??2773) 2 



(5.4) 



n: 



8p|p-fci||p + fc 2 | 

Xe - i |p-fcl|( r 11-»72) e -i|p+fc2|(l)2-'j3) e -ip(»73-!fl) _|_ . . . 

In the high momentum region of the loop momentum, 
the phase factor will be independent of the integrated 
momentum and the integral will be logarithmically di- 
vergent. As with the self-energy case before, these ar- 
guments can only demonstrate under what conditions a 
divergence can occur. Whether these logarithmic terms 
cancel among each other or whether the resulting diver- 
gence can be removed by a counterterm requires per- 
forming the full integration and summing all the rele- 
vant products of Wightman functions. However, since 
we have already seen that the self-energy graphs exhibit 
a pathological behavior in the a-vacuum, we shall not 
study these vertex corrections further here. 



VI. CONCLUSIONS 

The preceding discussion shows that a class of linear 
divergences from loops with two propagators — and loga- 
rithmic divergences from loops with three propagators — 
generically appears in any interacting theory in an a- 
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vacuum. These divergences arise from the form of the a- 
vacuum propagator, which is determined by the free field 
Hamiltonian, and not on the detailed form of the interac- 
tion. What the Schwinger-Keldysh formalism allows is a 
precise statement of the problem of an interacting theory 
in an a-vacuum. With a high-momentum cut off, we can 
find all the terms that diverge linearly with this cut off 
and analyze their dependence on the momenta external 
to the loop. The resulting expressions are not cancelled 
by a set of de Sitter invariant counterterms of the same 
form as those in the original Lagrangian. The appearance 
of a-dependent prefactors also shows why such terms do 
not plague the Euclidean vacuum. 

The fact that the divergence originates from high mo- 
mentum modes and the form of the a-dependent prefac- 
tors provides the basis of a simple heuristic explanation 
for the divergence. 1 The number density of Euclidean 
particles per unit volume at the time that the interac- 
tions are turned on is 



V-Uala^afla) 



N z 



1 



a? . From the perspective of the Euclidean 

k 



since a^(f] 

vacuum, the a-vacuum looks like a distribution whose 
occupation number is given by n|. Note that n2 is actu- 
ally independent of k. If we then replace the a-dependent 
terms in the propagators with the factors rip maintaining 
the momentum labels, then among the many divergent 
terms contributing to Ai occurs the expression, 



Ai = - 



r d jt r. 

Jva V J 



d 3 p 



A 2 v <WfMI ; 

kr) 64tt 3 T]l J VQ jf J p\p-k\ 
x [(„g + l)(n| + l)nZ^ k - + 1)] 

x sin[[p + k — \p— k\](rj — r/)l 

(6.2) 



in the massless conformally coupled case. This expres- 
sion resembles a 'gain minus loss' process in the a- 
background — for example, one part describes the cre- 
ation of two particles from one while the other describes 
the creation of one from the annihilation of two. Since 
rip is constant, nothing suppresses the large p divergence. 

This divergence is only present in a true a- vacuum and 
not in a 'truncated a-vacuum' — a state that is set equal 
to a Euclidean vacuum above some scale \k\ > M [TT| . 
For a truncated a-vacuum, the n^'s vanish above M so 
integrals such as Eq. (|6.2(l become finite. The largest con- 
tribution to the change in the number operator scales as 
X 2 M. These truncated a-vacua have no divergences, al- 
though exactly how the state evolves may depend on how 
the truncation is implemented. Using our formalism, it 



becomes possible to study how one of these vacua evolves 
during inflation and the amount by which it would alter 
the appearance of angular power spectrum of the cosmic 
microwave background radiation |l8| . 

To conclude then, interactions destabilize the a-vacua. 
We have calculated a physical quantity, the conformal 
time rate of change of the number of Euclidean mode 
particles in the a-vacuum, and found that it diverges. 
This divergence reflects a physical pathology of an inter- 
acting theory in a true a-vacuum. 
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APPENDIX A: REGULARIZATION OF 
LOGARITHMIC DIVERGENCES 

The self-energy diagrams in the Fig. [3] and Fig. ^con- 
tain two classes of divergences in the a-vacuum. One 
class diverges linearly in the magnitude of the spatial 
loop momentum. The dependence of this divergence on 
the external momentum flowing through the graph is not 
of the same form as that appearing from the insertion of 
one of the available counterterms. In this sense, these di- 
vergences can not be renormalized and indicate a patho- 
logical feature of an interacting theory in the a-vacuum. 
However, these linearly divergent terms vanish in the Eu- 
clidean limit. 

The second class of divergences exist in both the Eu- 
clidean and the general a-vacua. These divergences de- 
pend logarithmically on the conformal time and can, at 
least in the massless conformally coupled case where the 
calculation simplifies, be removed by a constant mass 
countcrtcrm. It is important to establish the renormal- 
ization of this type of divergence not only to show that 
the evolution of the number operator is finite for an in- 
teracting scalar field in the Euclidean vacuum, but also 
in the a-case. If we wish to consider a 'truncated a- 
vacuum', one which is cut off at some high energy scale 
such as the Planck mass 0] 1 we thereby remove the lin- 
ear divergence of Eq. Ij4.43|) . but the logarithmic diver- 
gence is still present in Ai and A®. In this appendix, we 
demonstrate how to renormalize this divergence. 

In the massless conformally coupled case, loop integral 
over a-vacuum propagators of Eq. H4.31(l yields 



1 We would like to thank Dan Boyanovsky for suggesting this ki- 
netic interpretation 



k.i 



2 



-ifc(?)l-772 ) 



a+a* gife(r;i— J72) 



1]1 - V2 
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+^-Kimm) V ~ e a } 

2 r)i+ r ]2 



ITT 



k 

ITT 

T 

2n 



e a N^rj 1 ri2 sin krji 



g 2a* e -ikri2 



,—ikrii _ e 2a 2ifc7ji"| 



— e a ' N^r]^ sin kr] 2 [e 

sinfc(?7i - rj 2 ) sin k(r]i + 7/ 2 ) 



+=j-(A-k)e a+a N^( mm ) 2 



(Al) 

m -m m + m 

The final term is the linearly divergent term. Among the 
remaining terms, only the first produces any logarithmic 
divergence when integrated over the conformal time. In 
the case of either self-energy contribution, Eq. (|4.3U|) or 
Eq. (14.32(1 . we integrate r\ 2 from tjq to r\\ {r\\ — > 77 for the 
Ai graph) and encounter a singularity from the (rjx — 
772) denominator of imaginary part of the first term in 
Eq. iTHll . 



(A2) 



2 v 



vi - m 

Inserting this result in the self-energy contributions yields 
A 2 V 



Ai 



krf 167T 



r(14 
drf 



)K\Uk(vo)Y 



— cos k(rj — tj') 



Vo iiv - v'j 
x\ (l + e a+a *) cos k(r]-T]') 



f e at e ik(r)+r]') _|_ £ a* e -»k(?j-fV)j 



A 2 y 
' fcr? 32tt 



(1 



??0 



— cos fe(r/ - rf) 
x|(l + e a+Q *)cosfc(?7-7;') 



+ ■ 



(A3) 



and 
-4o = 



A 2 



(1 



k 2 jf 16tt 2 

/"' d»7i Z"' 71 di]2 cos fe (771 - 772) 
x / / VaVb 

x sinfc(?7 Q - 771) (l + e a+a ) cos k(rj b - rj 2 ) 



_ e a e ik(Vb+V2) 



e e 



(A4) 



In both of these equations, the ellipses indicate terms 
which do not diverge logarithmically in the conformal 
time. 

Both Eq. (|A3|) and Eq. (|A4|) contain divergent integrals 
of the form 



di] 2 



(A5) 



where q is k, 0, or 
variable, 



Eq. <| A5|> becomes 
rvi e 2i ^ 2 dr] 2 



—k. Changing to a dimensionless 



Vi 



o^iqm r 1- ^ e~ 2iqVir dr 



Hi) 



V2(vi - m) 



Vi Jo 



(r - l)? 



(A6) 



(A7) 



We can regularize this integral by inserting a factor of ? ,e 
in the integrand and then extract the pole as e — > 0, 



' m e~ 2iqvir r e ~ 1 dr 1 
r — 1 e 



finite. 



(A8) 



Thus, 



— ^- = + finite. (A9) 

V2 (m - 772) Vi e 



'/o 



Applying this regularization to the self-energy graphs 
gives 



A, 



-(1 



A 2 V 1 

77 6 167T 2 e 

rEf„ \|2 



n+a 



Vo 



(A10) 



and 



2A 2 V 1 



if IQtt 2 £ 

4 Tm{us(va)ur(mM{vb)ur(rn)} 



+ ■ 



(All) 



In deriving Eq. I(A11(I we have used the fact that the 
operator D a b(v) is symmetric in 77a and 77b. Both poles 
are removed by the appropriate mass counterterm graphs 
Ci and C given in Eq. (|4.29|l and Eq. 14.33fl , respectively, 
when 



5m = 



1 A 2 



e 16tt 2 



a-\-a 



Ni 



(A12) 
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